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ABSTRACT 

Using multi-site spectroscopic data collected from three sites, the frequencies and 
pulsational modes of the 7 Doradus star HD 12901 were identified. A total of six 
frequencies in the range 1-2 d _1 were observed, their identifications supported by 
multiple line-profile measurement techniques and previously-published photometry. 
Five frequencies were of sufficient signal-to-noisc for mode identification and all five 
displayed similar three-bump standard deviation profiles which were fitted well with 
(7,m)=(l,l) modes. These fits had reduced % 2 values of less than 18. We propose that 
this star is an excellent candidate to test models of non-radially pulsating 7 Doradus 
stars as a result of the presence of multiple (1,1) modes. 

Key words: line: profiles, techniques: spectroscopic, HD12901, stars: variables: gen- 
eral, stars: oscillations 



1 INTRODUCTION 

The identification of the geometry of a 7 Doradus type pul- 
sation remains one of the more difficult spectroscopic fields. 
However, it is a powerful way to characterise the interior 
structure of a star and improve stellar pulsation models. The 
pulsations of 7 Doradus stars are particularly of interest as 
they propagate through the deep layers of a star. The flux- 
blocking mechanism at the base of the surface convective 
zone (Guzik et al. 2000) is held accountable for the origin of 
the pulsations which are g-modes, where the restoring force 
is gravity. 

Briefly 7 Doradus stars are slightly evolved late-A to 
late-F stars at the cooler end of the classical instability strip 
that display high-order g-mode pulsations with frequencies 
on the order of 1 cycle-per-day. For a full definition see Kaye 
et al. (1999) and recent reviews of the field are in Kaye 
(2007) and Pollard (2009). Currently there are less than 100 
bright bona fide 7 Doradus stars known (see list in Henry 
ct al. 2011) and a handful of 7 Doradus/J Scuti hybrid stars 
(Henry & Fekel 2005; Uytterhoeven et al. 2008), with a fur- 
ther 100 7 Doradus and 171 hybrid stars thus far reported 
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by the Kepler mission (Grigahcene et al. 2010; Uytterhoeven 
et al. 2011). 

Now is truly the age of satellite photometry and no 
ground-based methods can compete with the long unin- 
terrupted datasets that satellites such as C0R0T and Ke- 
pler produce. The Fourier frequency spectra of such studies 
show the high-levels of precision only obtainable from space 
(e.g. at least 840 frequencies found in HD 49434 using the 
C0R0T satellite; Chapellier et al. 2011). However, to further 
our understanding of g-mode pulsations we need more data 
than frequencies alone. Several successful techniques using 
ground-based multi-colour photometry can be employed to 
determine the number of nodal lines (I) in a star. The full 
mode-identification of a star, that is finding the I and also 
the number of nodal lines passing through the pole of a star 
m, remains the sole domain of spectroscopy. It is hoped that 
modelling of the frequencies and modes of a star will allow 
determination of the n, number of shells interior to the stel- 
lar surface. 

The spectroscopic study of 7 Doradus stars relies on the 
collection of a large amount of high-resolution and preferably 
multi-site data of sufficient signal-to-noise for classification 
of the pulsations. Improvement is gained by the use of data 
from multiple sites to reduce daily aliasing patterns in the 
Fourier spectra. Such observational campaigns take many 
months to years in order to collect a sufficient number of 
spectra. 

When such datasets are compiled for 7 Doradus stars, 
the pulsational frequencies and modes can be examined and 



Table 1. Observation log of the multi-site data for HD 12901 

Observatory Telescope Spectrograph 



Observation range AT (d) # spectra 



MJUO 1.0m HERCULES Feb 2009-Scp 2011 935 478 

La Silla 1.2m CORALIE Nov 1998-Nov 2002 1166 53 

McDonald 2.1m SES Sep 2009-Oct 2009 8 60 



compared to those from photometry. There are still only a 
handful of 7 Doradus stars with full mode identifications 
and our immediate goal is to classify as many as we can. 
The results of these mode identifications can be used to im- 
prove pulsational models (such as those of Grigahcene et al. 
2012; Townsend 2003) by providing information about the 
amplitudes of excited modes and also to start to identify 
patterns within the class. 

This paper focuses on the 7 Doradus star HD 12901 
outlining the observations made and reduction procedure 
in Section 2. The spectroscopic frequency analysis for the 
individual observatories and the combined data results are 
described in Section 3 with each method tested. Section 4 
describes the reanalysis of white-light and seven-colour pho- 
tometry taken of this star. The mode identification of the 
five identified frequencies follows in Section 5. Finally Sec- 
tion 6 discusses the findings and their implications for future 
work in this field. 



2 OBSERVATIONS AND DATA TREATMENT 

The major findings of this research are the identification of 
the frequencies and modes using spectra obtained at three 
observing sites (Mt John University Observatory (MJUO), 
La Silla Observatory and McDonald Observatory) and are 
summarised in Table 1. 

Spectra were reduced according to the standard local 
spectrograph software then they were normalised, contin- 
uum fitted and order-merged by the authors using a semi- 
automated MATLAB routine. The full spectra were cross- 
correlated for each site using a scaled delta function routine 
(Wright 2008, Wright et al. 2007). This produces the line 
profiles used for pulsational analysis. The signal-to-noise of 
a cross-correlated line profile is much higher than a single 
spectral line and thus pulsations are easier to extract. The 
line profiles for each site were combined (scaled, shifted in 
velocity space and continuum fitted) to produce a single 
dataset of consistent line profiles. This was done post-cross 
correlation to maximise the number of available lines used 
for each spectrum, as different spectrographs operate over 
different wavelength ranges. 

The line profiles were analysed in FAMIAS, a pulsational 
frequency mode-identification toolbox (Zima 2008), using 
techniques applied from Zima et al. (2006). The Fourier 
spectra for the Pixel-by-Pixel (PbP) technique and that of 
the 0th-3rd moments (Balona 1986; Briquet & Aerts 2003) 
were used to identify frequencies. These frequencies were 
then analysed to determine mode identifications using the 
Fourier Parameter Fit method (Zima 2009). 

Additionally, the software package SIGSPEC (Reegen 
2007) was used to compare the frequency selection method 
of FAMIAS. SIGSPEC performs a Fourier analysis of a two- 
dimensional dataset and selects frequencies based on their 
spectral significance. The frequencies obtained with SIGSPEC 



can be regarded as more precise than that of the highest- 
peak direct selection as, after producing the Fourier fre- 
quency spectrum, the peak with the highest spectral sig- 
nificance is selected. Spectral significance includes the anal- 
ysis of the false-alarm probability to remove frequency peaks 
caused by irregular data sampling or noise in the data. 

Both FAMIAS and SIGSPEC were used to re-analyse pho- 
tometric data, originally published in Handler (1999), Eyer 
& Aerts (2000) and Aerts et al. (2004), as further insight 
into the frequencies of the pulsations. To identify the I val- 
ues of the modes, the frequency amplitude ratio and phase 
differences method, based on pre-computed grids of models, 
were used (Balona & Stobie 1979; Watson 1988; Cugier et al. 
1994; Daszynska-Daszkiewicz et al. 2002). This was done us- 
ing the photometric analysis toolbox in FAMIAS (Zima 2008). 



3 SPECTROSCOPIC FREQUENCY ANALYSIS 

In total 591 spectra from the three sites were of sufficient 
quality to be analysed. The data spanned a total of 4667 
days, just over 12 years. Each dataset was analysed inde- 
pendently for frequencies and they were then combined. By 
doing this analysis we can comment on the limitations of 
single-site data and the extent to which multi-site data re- 
duces aliasing. 



3.1 Frequency Identification For Each Dataset 

3.1.1 MJUO 

This is the largest single-site dataset and was analysed most 
extensively. First the cross-correlated line profiles were anal- 
ysed in FAMIAS. The frequency list is in Table 2. The PbP 
technique found eight possible frequency peaks including 
a one-cycle-per-day frequency (/ m5 ), likely an alias. The 
zeroth moment was too noisy to extract any frequencies. 
The first and third moment Fourier peaks each only showed 
one frequency clearly, 1.186 d _1 and 1.184 d _1 respectively, 
which is the same frequency given a conservative uncertainty 
estimate of ±0.001 (uncertainties are dealt with more for- 
mally in Section 3.3). The second moment yielded 0.144 d _1 
as the only viable frequency. Though many more frequencies 
may be evident in the data, further extraction of frequencies 
above the noise was difficult with any measure of certainty. 

The same moments were then tested in SIGSPEC for fur- 
ther analysis. SIGSPEC found 20 frequencies above a spec- 
tral significance level of 5. It is likely that noise in the data 
is causing many misidentifications rather than the possibil- 
ity that these could all be real, given the noise level in the 
Fourier spectra. The frequencies with spectral significance 
greater than 15 are reported for moments one to three in 
Table 2. From the table it is clear that f m 2 to f m s (excluding 
/ms) are viable frequencies found using multiple methods. 
The addition of further datasets from other observation sites 



should improve the signal-to-noise and reduce any aliasing, 
particularly 1-day aliasing that occurs in the data. 

3.1.2 La Silla 

The spectra taken on the CORALIE spectrograph have pre- 
viously been analysed and published by Aerts et al. (2004) 
and De Cat et al. (2006). Frequencies of 1.04±0.28 d _1 and 
1.30 ± 0.30 d _1 were reported. 

Although the small size of the dataset made it difficult 
to analyse, frequency peaks were identified at 0.005 d" 1 and 
1.26 d" 1 in the PbP measure, 0.997 d" 1 and 3.01 d _1 in 
the first moment and 0.997 d" 1 and 1.999 d" 1 in the third 
moment. Peaks that occur very close to integer values are 
less reliable as they are likely 1-day alias patterns from the 
observation times artificially amplified by data sampling ef- 
fects in the window function. No identifiable peaks occurred 
in the zeroth and second moments. Only the 1.26 d _1 can 
be identified with any confidence. It is notable that a peak 
was apparent in all the above Fourier spectra at 0.23 d _1 . 

SIGSPEC was better able to distinguish the 1-day aliases 
and identified frequencies of 0.26 d -1 in the zeroth moment 
and 0.23 d _1 in the first and third moments. These had sig- 
nificances of 3.3, 3.5 and 2.7 respectively, which are usually 
regarded as too small to be significant. 

3.1.3 McDonald 

The McDonald data set comprised of 60 observations taken 
in September and October 2009. The data were considered 
to be particularly useful as having time intervals that overlap 
some of the MJUO observations, providing some independent 
confirmation of the line profile variation. As a stand-alone 
data set the shorter wavelength range of the spectrograph, 
and the larger regions of telluric lines, reduced the number of 
stellar lines that could be cross-correlated from several thou- 
sand to around 100 lines. This has a significant impact on 
the signal-to-noise of the final line profile. This, and the low 
number of observations, meant frequencies detected in this 
dataset alone were unconvincing. Despite this, it is noted 
that frequencies near 1.38 d" 1 and 1.24 d _1 are present in 
the PbP analysis and a peak near 1.67 d _1 is visible in the 
first moment Fourier spectrum. 



3.2 Frequency Identification Of Combined Data 
For Each Method 

The cross-correlated line profiles of the three datasets 
were combined post cross-correlation to form a single 
dataset. The line profiles are displayed in Figure 1. Com- 
bining the line profiles instead of the individual spectra was 
the best way to preserve the signal-to-noise in the MJUO and 
La Silla datasets which had many more lines than the Mc- 
Donald spectra. There was a small shift (6 kms -1 ) in the 
systemic velocity between the McDonald and the other two 
datasets. The McDonald spectra were shifted to ensure con- 
sistent data. The line profiles were scaled to have the same 
equivalent widths, a parameter shown to remain constant 
by the lack of variation in the zeroth moment in the MJUO 
data. The Fourier spectra of the MJUO data alone and the 
combined datasets are compared in Figure 2. 
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Figure 1. Scaled line profiles of the three datasets. McDonald 
(red), La Silla (blue) and MJUO (black) observations are shown 
with the mean line profile (white). All datasets scale well to show 
a consistent line profile. The data are provided as an electronic 
file (see Supporting Information). 



The spectral windows for the MJUO data alone and all 
three data sets are very similar, but, as expected, the ad- 
dition of the multi-site data reduces the amplitude of the 
secondary peaks of the identified frequency. From Figure 
2 it is clear that the addition of the extra data increases 
the amplitude of the Fourier peaks. This is an effect of not 
only the reinforced signal found in the other sites spectra, 
but also an artificial amplification due to the increased base 
noise level. The frequency peaks remain in nearly the same 
positions and show approximately the same repeating pat- 
tern. The one-day aliasing, although still present, has been 
reduced to make the true peaks more evident. In general, 
the frequency identification becomes clearer with the addi- 
tion of the multi-site data, despite it being from small sets 
compared to the MJUO data. 

3.2.1 Pixel-by- Pixel (PbP) 

The line profiles showed many periodic variations in each 
pixel. The ten frequencies found in the combined dataset 
are listed in Table 3 and the Fourier spectra displayed in 
Figure 2. The first two frequencies found were the fa and fz 
from Aerts et al. (2004) from photometry and their /i shows 
up as our f P 4. Further frequencies / P 3 and / p s have credible 
three bump variations in their line profiles. It is likely that 
f P 6 is a one-day alias frequency arising from the dominance 
of one single site (MJUO contributes about 80% of the ob- 
servations). The frequencies f p r, f p $ and f p io showed clear 
three bump standard deviation profiles each with smooth 
phase changes. This increases the likelihood that they are 
real frequencies, although the possibility remains that they 
are from residual power from a previous poorly-defined fre- 
quency. Finally, f p g showed a clear four bump standard de- 
viation profile at 3.14 d~ and 4.14 d _1 frequencies when 
tested. As none of the previous frequencies showed a four- 
bump pattern in the variation it is more likely this is a real 
frequency rather than a misidentification of an earlier fre- 
quency. 

A check on the independence of the frequencies was 
done by prewhitening the spectra by 1-day aliases and 1- 
/ aliases of the first five identified frequencies. Only f p s ± 1 
had clear standard deviation profile variations resembling 
pulsation and none of the alternate identifications altered 



Table 2. Frequencies of HD 12901 found in the MJUO data. Pixel-by-Pixel (PbP) frequencies found using FAMIAS and moments and 
significances (Sig) found using SIGSPEC are listed. The amplitudes of the PbP frequencies are scaled to the first identified amplitude. 
Note some frequencies have been identified as 1 — / aliases of the PbP frequency. The frequency f m s shows identifications near 1 d _1 , 
likely from the data sampling. 
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21 
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23 


fm3 


1.681 
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Figure 2. Effect of the multi-site datasets on the Pixel-by-Pixel Fourier spectra, the MJUO data on the left and the combined data on 
the right. Each series shows the window function for the data at the top then successive pre-whitening of frequencies f v \ to f p §. The 
Fourier spectra for the Pixel-by-Pixel analysis of the combined data continues to f p io- The frequencies are listed in Table 3. 



Table 3. Frequencies from the Pixcl-by-Pixel analysis of the com- 
bined dataset. The uncertainty estimate for the frequencies is 
±0.0002d~ 1 as described in Section 3.3. Frequencies with a strike- 
though were discarded as described in the text. 



ID 


Freq. 


Period 


Amp. 


Phase 


Variation 




(d- 1 ) 


(d) 


(rel. to / p i) 




Explained 


fpi 


1.3959 


0.7164 


1 


0.286 


13% 


fpi 


1.1863 


0.8430 


0.72 


0.205 


21% 


fp3 


1.6812 


0.5948 


0.93 


0.267 


30% 


fpi 


1.2157 


0.8226 


0.93 


0.265 


39% 


fp5 


1.5596 


0.6412 


0.80 


0.230 


45% 


fp6 


1 0004 


0.9996 


0.55 


0.158 


48% 


fp7 


1.2465 


0.8022 


0.50 


0.143 


50% 


fp8 


1.2743 


0.7848 


0.53 


0.152 


53% 


fp9 


3.1302 


0.3186 


0.36 


0.102 


54% 


fpW 


2.5803 


0.3876 


0.37 


0.106 


55% 



the identification of following frequencies, indicating that 
they are all independent. 

A second check of the variability of the line profile with 
different frequencies involves phasing the data to the pro- 
posed frequency and examining the structure. To see the 
changes most clearly this was done using the residuals of 
the line profiles after subtraction of the mean. The frequen- 
cies / p i to /p 5 in Figure 3 the residuals on the left side of 
each plot show the 'braided rope' structure typical of non- 
radial pulsation, but for the other frequencies this was less 
clear. If we require pulsations to produce regular changes 
in both red and blue halves of the line profile then we must 
discard f P 6, fpi and f p g as candidate frequencies leaving just 
/p8 and /pio to investigate further. 

An investigation into the phase coverage of the data for 
each frequency can also be used as an indicator of the relia- 
bility of the frequency. All of the frequencies f p i — /pio were 
well covered in phase space except / p6 , which was poorly 
sampled between phases 0.2-0.3 and 0.5-0.8 due to the dom- 
inance of a single site in the data. 



cies found in the highest peak method then reproduced PbP 
frequencies f P 3 and f p j as //g and //io- 

The Fourier spectra of the second moment also showed 
a few promising frequencies once the first peak at 0.0006 d _1 
was removed. The frequencies f S 2 to f„4 match to, or match 
to one-day aliases of, frequencies found in the PbP dataset. 
It is likely that f a6 is a misidentification of f p2 . The second 
section of Table 4 catalogues the frequencies found. Both 
the highest peak and PbP frequencies accounted for around 
the same amount of variation as they identified nearly all 
the same frequencies. 

Frequency peaks found in the third moment are pre- 
sented in the final section of Table 4. The first few frequen- 
cies found using this method were generally similar to those 
found in the PbP method except the double identification 
of ft2- This is usually the result of small errors in the iden- 
tification of a strong frequency, leaving residual power in 
the Fourier spectrum. It is likely that this then impacts the 
following sequence fts to /no which all seem close frequen- 
cies to other identifications, and possibly one-day aliases. 
The PbP frequencies identified explained more variation, ac- 
counting for a total of 83% of the line profile variation. 

Overall the PbP variations have been shown to remove 
the variation from the data more accurately, demonstrated 
by the higher percentage of variation removed for each mo- 
ment. This leads us to conclude that the frequencies found 
using the highest peaks in the moment methods are less re- 
liable. It is notable that the SIGSPEC identified frequencies 
are more reliable than the highest peak selection. 

As an additional test of the robustness of the derived 
frequencies, a synthetic dataset, modelled on the time spac- 
ings and modes identified in Section 5, was created to com- 
pare with the PbP derived frequencies. The synthetic line 
profiles were phased and plotted to compare with the real 
data in Figure 3. The synthetic profiles strongly match the 
observed profiles, strengthening our confidence in the fre- 
quency identification. 



3.2.2 Moment Analysis 

The zeroth moment analysis shows few peaks that are simi- 
lar to those from other methods. There is a cyclic variation 
around 1-cycle per day but with the 1.0 d" 1 peak missing 
due to the combination of multi-site data. The small ampli- 
tude of variations in the second moment indicates we have 
only small periodic temperature variations to account for in 
the line profiles. 

The first moment Fourier spectra showed clear peak fre- 
quencies similar to the PbP method. These are shown in the 
first section of Table 4. The amount of variation of the first 
moment was moderately well explained by the selection of 
the highest six peaks in the Fourier spectrum. Beyond ff^ 
the improvements to the explanation of variation are too 
small to be conclusive. When the PbP frequencies were ex- 
tracted from the first moment data (in approximately the 
same order as the highest peak frequencies) a better fit to 
the variation was found. With the six frequencies, f p \ — / p5 
and f p s, 76% of the variation was removed. The frequencies 
found using SIGSPEC closely matched the first four frequen- 



3.3 Frequency Results 

The frequencies chosen to proceed with in the analysis were 
those originally identified in the PbP results, f p \ to / p5 and 
/p8 as noted in Table 3. This choice was made based on the 
prevalence of these frequencies through the other techniques 
used, the decreased sensitivity of the PbP technique to 
asymmetric variations and generally observed higher signal- 
to-noise of the PbP frequencies. 

It is clear from the above sections that, due to the inclu- 
sion of multi-site data, 1-day aliasing was not a big problem 
when identifying frequencies present in these data. 

Identifying the uncertainty range on the frequencies 
found is a complex task. When using SIGSPEC we can use 
the relation of Kallinger et al. (2008): 

*(/)= )=r= v « 

T * ^/sig(a) 

where T is the time base of observations in days and sig(a) is 
the spectral significance of the frequency. For the frequencies 
/pi to /pio identified in the first moment method this gives 
an uncertainty estimate of ±0.0002. Additionally we used 
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Figure 3. Line profile residuals phased on the frequencies iden- 
tified /pi (top) to / p 5 (bottom) compared with synthetic pro- 
file TU„ ^ r .f+ nk™^. ~Ur,™r„A P 1^K+ Ir. 



Table 4. Resulting frequencies from the First, Second and 
Third Moments found using peaks in the Fourier spectra and 
in SIGSPEC 



First Moment 





Highest Peak 


SIGSPEC 


ID 


Freq. 


Var. 


Freq. 


Sig. 


//i 


1.1835 


16% 


1.18624 


26 


ff2 


1.6100 


32% 






ff3 


1.4017 


42% 


1.39591 


26 


ft* 


1.2247 


52% 


1.2156 


28 


1.2154 




2.2273 


21 


ffs 


0.6490 


57% 




22 


ff6 


1.5627 


63% 


1.5596 


20 


fir 


0.7559 


66% 






ffs 


1.5179 


69% 






ff9 






1.6812 


22 


fno 






0.2436 


22 


//ii 






2.1362 


17 


//12 






0.0002 


17 


/ /13 






2.7789 


15 


ffli 






0.6736 


15 




Second Moment 






Highest Peak 


SIGSPEC 


ID 


Freq. 


Var. 


Freq. 


Sig. 


fsl 


0.0006 


7% 






fs2 


0.2159 


22% 






fa3 


0.3767 


38% 








1.3962 




1.3959 


19 


fsi 


1.5596 


46% 


1.5596 


20 


fs6 


2.1786 


57% 






fs7 






0.1622 


25 


fs8 






0.9449 


23 


fs9 






0.9061 


19 


/slO 






0.7382 


19 


/all 






2.6809 


17 




Third Moment 






Highest Peak 


SIGSPEC 


ID 


Freq. 


Var. 


Freq. 


Sig. 


fa 


1.3987 


17% 


1.3956 


31 


ft 2 


1.1864 


40% 


1.1862 


25 




2.1833 








ftz 


1.6810 


49% 


1.6809 


23 


fH 


1.3489 


58% 






fts 


1.1661 


64% 






fte 


2.5624 


32% 


0.5597 


21 


ftr 


2.7952 


41% 






fts 


1.1716 


49% 






ft'.) 


3.5169 


58% 






/tio 


1.2658 


64% 






/til 






1.2156 


27 


/il2 






2.1923 


20 


/tl3 






0.2437 


18 



the estimation from Montgomery & Odonoghue (1999) who 
propose 



6 1 a(m) 



(2) 



NnT A 

Here N is the number of observations. The value a(f) is 



Table 5. Frequencies found in the HIPPARCOS white-light pho- 
tometry. The first three frequencies match well with f p 4, f p s and 
fp7- 



Frequency (d 1 ) 


Uncertainty 


Significance 


1.2158 


±0.0003 


7 


1.5525 


±0.0003 


6 


2.2449 


±0.0003 


6 


2.5665 


±0.0004 


5 



a one-sigma uncertainty with an amplitude (A) root-mean- 
square deviation a(m). This method was used to derive (an 
underestimate) of the uncertainties of the pixel-by-pixel fre- 
quencies. The results for the frequencies f p \ to f p5 ranged 
from ±0.00013 to ±0.00018. Given the observed differences 
between the frequencies identified in the PbP and the mo- 
ment methods we view ±0.0002 to be a good estimate of the 
uncertainty. 

The data were also analysed for possible frequency 
combinations. It was found that /ps(1.559) ± / P 3(1.681) = 
3.240d _1 which is possibly the same as / P 7(1.246) and 
strengthens our earlier removal of this frequency. 

It is clear from the above sections that many frequen- 
cies extracted from the data are robust as they appear in 
multiple methods with clear line profile variations. The first 
five frequencies in the PbP analysis appear in almost every 
analysis method. The high signal-to-noise in the line profiles 
for these frequencies makes them suitable for mode identi- 
fication. The frequency f p s, although showing evidence for 
being an independent frequency, showed a misshapen stan- 
dard deviation profile, which meant there was not enough 
certainty in this frequency identification to proceed with a 
mode identification. The final PbP frequency f p io showed 
no evidence of appearing in other analysis methods and was 
rejected for mode identification. 

Even with all these frequencies identified, variations in 
the data remain. The PbP method only removed 50% of 
the variation and the moment methods up to 80%. Some of 
the remaining variation is due to the noise in the data and 
possible slight misidentifications of the frequencies, but it 
is likely that there remain multiple unidentified frequencies 
below the detection threshold of the present data. The pres- 
ence of these further frequencies is expected for 7 Doradus 
stars as they have dozens and sometimes even hundreds of 
frequencies identifiable in photometry. The large number of 
similar amplitude frequencies make this star challenging to 
study. However already we are able to see more frequencies 
than previously identified in spectroscopy, which makes this 
a promising star for further study. 



4 PHOTOMETRIC FREQUENCY AND MODE 
IDENTIFICATION 

Two photometric datasets were analysed for frequencies. 
The first was a time series of white light observations from 
the satellite HIPPARCOS (Perryman & ESA 1997). Specifi- 
cally these were taken in the H p filter. The data span a pe- 
riod of 1166 days from November 1989 to February 1993 dur- 
ing which 122 measurements were taken. The second dataset 
was Geneva photometry taken on the 0.7m Swiss telescope 



at La Silla over a period of 25 years from 1973 to 1997. A 
total of 174 observations for each filter in this time. This 
dataset has been extensively analysed in Aerts et al. (2004). 
We present our re-analysis to complement the spectroscopic 
results above. 

4.1 HIPPARCOS Photometric Frequencies 

The Fourier analysis of the HIPPARCOS data measurements 
have a range of 0.09 H p magnitudes with an average uncer- 
tainty of 0.1 Hp magnitudes. It was difficult to identify clear 
frequencies from the Fourier spectra, but the highest peak 
was at 1.2701 d~ . Using SIGSPEC we were able to identify 
the frequencies and their significances as shown in Table 5. 

4.2 Geneva Photometric Frequencies 

The seven filters of the Geneva photometry provide us with 
enough information to extract frequencies and identify the 
/ modes. Each of the seven filters showed very similar vari- 
ation and this was reflected in the individual Fourier spec- 
tra. Frequency peaks were observed at 1.218 d" 1 , 1.396 d" 1 
and 0.0907 d" 1 (or one-day aliases) in most filters. To for- 
malise the frequency identification we used SIGSPEC to iden- 
tify the frequencies and their significances as shown in Table 
6. These frequencies match well to those previously found in 
the same data as is discussed further in Section 6. 

4.3 Geneva Photometry Modes 

The identification of the I value of the spectroscopic fre- 
quencies / p i to /p5 was attempted using the seven filters in 
the Geneva photometric data. The identification was done 
using the amplitude ratio method in FAMIAS with modes 
from I = 1 to I = 3 tested in all filters. All five frequencies 
were found to be solely consistent with I = 1 modes distin- 
guished primarily by the first amplitude ratio. The results 
are largely identical to Aerts et al. (2004), Their Figure 6 
shows the unique fit to the I = 1 mode. This method, al- 
though not as powerful as spectroscopic mode identification, 
gives us independent support to the / = 1 mode fits of the 
star. 



5 MODE IDENTIFICATION OF 

SPECTROSCOPIC LINE PROFILES 



The mode identification was performed using the PbP 
frequencies f p i to / p s, as chosen in Section 3.3 (/ P 8 not being 
considered due to the distorted line profile variations). The 
modes of the individual frequencies were first identified and 
then a best fit, including all five frequencies, was computed. 
Initially the zero-point profile was fitted to determine the 
basic line parameters. The best fit is shown in Figure 5(a). 
The vsim was found to converge at 63.9 kms -1 with a \ 2 of 
103. Figure 4 shows the 95% confidence limit on the determi- 
nation of vsini which gives a range of 63.9 ±0.5 kms -1 . The 
zero-point fit set the initial line and stellar parameters of the 
space searched. The details of the parameters are given in 
Table 8 and the resulting best fit modes are given in Table 9. 



Table 6. Frequencies (d 1 ) found in each of the Geneva photometry filters and significances using SIGSPEC. 
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Figure 4. The 95% confidence limit (solid line) of the zero-point 
fit to the vsini. The dashed lines show the limits at the confidence 
level and the dotted line indicates the minimum and hence best 
fit value. 



Table 7. Rotational parameters and critical limits for a pulsating 
star with mass = 1.5 Mq, radius = 1.7 Rq, vsini = 64 kms — 1 
observed at an inclination of i = 30°. 



v ro t 

Trot 
frot 
v crit 

vsinicrit 

Lev it 



128 kms" 1 
0.67 d 
1.49 d" 1 
410 kms -1 
205 kms" 1 
9° 



star is not approaching critical rotational velocity. This is 
the case for all tested values of i in the range i = 10° — 90°. 
The rotational frequency of the star is dependent on the 
inclination, ranging from 4.28 d" 1 at i = 10° to 0.74d _1 at 
i = 90°. The rotational frequency for a 7 Doradus star is 
expected to be on the same order as the pulsation frequency 
and these values fall within this range. 



The unusually high values for mass and radius are discussed 
in Section 5.1. 

The best fit models found for each individual frequency 
are plotted in Figures 5(b) - 5(f). The results in Table 9 
and Figures 6(a) to 6(e) show that the mode identifications 
are well determined as there is no ambiguity in choosing the 
model with the lowest \ 2 - The inclination appears to lie in 
the region 20° — 45°, with the simultaneous fit giving a value 
of 27° for a \ 2 of 16.59. Formally the 95% confidence limit of 
this parameter gives 27° ^ili so ^ is poorly constrained by 
the mode identification. Given the vsini of the star, a value 
near 30° is plausible as described further in Section 5.1. The 
vsini value determined varies slightly between the fits as it 
is modified to change the width of the standard deviation 
profile. The measured vsini of the star is best determined 
from the fit to the zero-point profile given above. 

The detection of five independent frequencies and the 
occurrence of a large number of (1,1) modes makes this 
star an excellent candidate for further asteroseismic analysis. 
The prevalence of the (1,1) modes may indicate sequencing 
of the n-values, or the number of interior shells, ft also is 
possible that the (1,1) modes may be linked to the rapid 
rotation of the star. This is discussed in the context of all 7 
Doradus stars in Section 6. 

5.1 Rotation and Pulsation Parameters 

The results of a preliminary mode identification were tested 
using the rotation and pulsation parameter tools in FAMIAS. 
Using mass = 1.5 Mq and radius = 1.7 Rq (typical for a 
7 Doradus star), vsini = 64 kms -1 and various inclinations 
indicates the rotational parameters of the star. Shown in 
Table 7 are the results for i = 30°, which indicates that this 



The pulsation parameter tool can be used to get an indi- 
cation of the horizontal-to-vertical pulsation amplitude pa- 
rameter (k) and whether the pulsation and rotation frequen- 
cies fall within the mode determination limits of FAMIAS. 
The k value falls between the values 50-800 for frequen- 
cies /1 — fi, which is above 1.0 as expected for 7 Doradus 
stars. These values then show the natural pulsational fre- 
quencies to lie between 0.2 d _1 and 0.7 d _1 . Taking the 
ratio frot/ j co-rot gives values in the range 1-5 for the de- 
termined frequencies and modes (all m = 1). FAMIAS has 
an operational range of frot/ f co-rot ^ 0.5, being designed 
to deal with p-mode pulsations with much larger vertical 
pulsation components. The very high values of « for the g- 
mode pulsations in 7 Doradus stars and the high rotation 
rate of this star mean we are operating beyond this limit. To 
obtain a reasonable mode identification the mass and radius 
were extended to non-physical values for a 7 Doradus star 
to provide observable amplitudes of pulsation (see Table 8). 
However the mass and radius are solely used for the calcu- 
lation of the ratio k, and do not affect other aspects of the 
fit. Making these changes allowed us to consider pulsations 
with frot/ f co-rot values of 0.3-0.5 for the multi-frequency fit 
in Table 9. 

Including more physics to more accurately describe the 
effects of rotation of g-modes in current pulsation models 
would ultimately solve such problems. Townsend (2003) in- 
vestigated the effect of higher rotation on g-modes with high 
n and low I by increasing the Coriolis forces. It was seen that 
the identification of prograde modes (defined here as m > 0, 
in Townsend 2003 as m < 0), have smaller differences in the 
rotating scenarios than retrograde modes. Smaller values of 
m are also less affected. There is also a discussion of this 
in Wright et al. (2011) who estimate the true limits of the 



Table 8. Stellar parameters used in the mode identification. 
Fixed values for T e fj, [M/H] and logg were taken from Bruntt 
ct al. (2008) and vsini values from the ranges in Bruntt et al. 
(2008); De Cat ct al. (2006); Duprct ct al. (2005); Acrts et al. 
(2004); Eyer & Aerts (2000). 



Parameter 


Fixed 
Value 


Min 


Max 


Step 


Radius (solar units) 




1 


10 


0.1 


Mass (solar units) 




0.5 


50 


0.01 


Temperature (K) 


7193 








Metallicity [M/H] 


0.13 








log 9 


4.18 








Inclination (°) 







90 


1 


vsini (kms -1 ) 




60 


75 


0.1 



Table 9. Results of mode identification of all four frequencies 
individually after a least-squares fit is applied (lsf), and all four 
frequencies simultaneously (sim). 
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1.07 


0.81 
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2.90 


0.18 


fi sim 


(1,1) 






1.07 


0.71 


fa sim 


(1,1) 






2.17 


0.60 



frot/f co-rot ratio in FAMIAS for various in for a 7 Doradus 
range frequency and found FAMIAS could identify modes up 
to frot/f co-rot = 1 for a m = 1 mode. These results indicate 
that the mode identification is unlikely to be affected by the 
rotation of this star. We can be confident that using non- 
physical values for mass and radius affects only the observed 
amplitude of the pulsation and not the mode, itself. 

A convincing mode identification of the frequencies /1 — 
fa was obtained. 



6 DISCUSSION 

This study allowed a direct comparison between single-site 
and multi-site data. From this we can judge the usefulness of 
large single-site datasets. The findings show that although 
the addition of further sites increased the amplitudes of the 
frequencies, it also elevated the base noise level. Addition- 
ally it is clear from the window function spectra in Figure 2 
that the 1-day aliasing pattern is reduced but not entirely 
eliminated. The above leads us to conclude that the addi- 
tion of multi-site data is useful, but not required to extract 
frequencies from sufficiently large datasets. We must also 
require that the data from any one site have a sufficient 
number of observations to produce a balanced (well sam- 
pled in phase space) mean line profile in order to combine 
the cross-correlated line profiles with the highest precision. 

The frequencies found in the PbP method were the most 
reliable and were consistent with all the other analysis meth- 
ods. Past papers have analysed the photometric data from 



HIPPARCOS and also some multi-coloured photometry. In 
Handler (1999), the author finds a frequencies of 2.18 d _1 
and 1.2155 d" 1 (f p4 ). In this paper we found frequencies 
/ P 4, / P 5 and (2x) f p8 in the HIPPARCOS photometry. 

The multi-colour photometry was previously published 
by Aerts et al. (2004), who found frequencies equivalent to 
f P 4, fpi and a 1-day alias of f p s. In the same data we find f V 4, 
f p i and / p5 . The mode identification showed the frequencies 
f P 4, f P i and / P 3 to best fit Z = 1 modes, the same as found 
in this analysis. The same result was found by Dupret et al. 
(2005) who showed that the three frequencies from Aerts 
et al. (2004) are I — 1 modes when using time-dependent 
convection models. 

Aerts et al. (2004) also describe the spectroscopic 
dataset taken with CORALIE that was also used in this anal- 
ysis. The authors did not find any frequencies from the spec- 
troscopic data alone (it is, as the authors note, too small a 
dataset for spectroscopic analysis). The imposition of the 
fp4, fpi frequencies found in their photometric analysis did 
provide some harmonic fits with low amplitudes. 

A study applying the frequency ratio method to the 
frequencies found in Aerts et al. (2004) was done by Moya 
ct al. (2005), who found three models consistent with I = 1 
modes. The models have a T= 6760if, logp = 3.88 - 4.12 
and stellar ages around 2 — 3 Myr. 

Considering all the prior studies, all of the candidate 
frequencies we confirm (f p i to / P 5 and f p a) are well sup- 
ported. The mode identification is also partially confirmed 
by photometry. There is no suggestion of different frequen- 
cies observed in photometry and spectroscopy as there are 
from some 7 Doradus stars (see Maisonneuve et al. 2011; 
Uytterhoeven et al. 2008). 

A prevalence of (1,1) modes in this star, and in 7 Do- 
radus stars in general, is beginning to emerge. This includes 
two modes in HD 135825 (Brunsden et al. 2012), two in 
7 Doradus (Balona et al. 1996; Dupret et al. 2005), one 
in HD 40745 (Maisonneuve et al. 2011), one in HR8799 
(Wright et al. 2011), two in HD 189631 (Davie 2011) and 
two in HD 65526 (Greenwood 2012). This is possibly due 
to the large surface area covered in each segment of a (1,1) 
mode, meaning pulsations have larger amplitudes and large 
changes in amplitude across the stellar surface. This indi- 
cates an observational selection effect. Additionally Balona 
et al. (2011) discuss the prevalence of m = ±1 modes in a 
sample of 7 Doradus candidates from Kepler photometry. 
These authors assume the dominant frequency to be the 
rotational frequency and they then show further frequen- 
cies to be close to this value. This constrains the light max- 
ima to once per rotation cycle, requiring a m = ±1 mode. 
The occurrence of five of these such modes in this star sug- 
gests some physical linking between them. The identifica- 
tion of several (1,1) modes in this star led to an investiga- 
tion into the period spacings of the six identified frequencies 
(Table 3). The asymptotics of oscillation theory (Tassoul 
1980) predicts a characteristic period spacing for high-order 
g-modes of the same low degree (I) for sequential values of n. 
An investigation into the period spacings of the PbP identi- 
fied frequencies shows that the spacing between f p i-f P 2 and 
fpi-f P 3 to be close (0.1266d and 0.1215d respectively). This 
suggests they could be subsequent values of n if we allow 
our frequencies to vary by ±0.003 d _1 . The frequencies f P 4, 
f P 5 and fpg however do not fit with this spacing and addi- 
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(f) Mode identification of f p5 = 1.5596 d" 1 (x 2 = 5.5). 



Figure 5. The fit (dashed) of the mode identification to the mean line profile, variation and phase (solid) of the five identified frequencies. 
All have been identified as (1,1) modes with the labelled x 2 for the best fit. An indication of the maximum uncertainty is given on each 
plot. 



tionally it may be that the close spacing of f P 2 and f P 4 is 
inconsistent with current theoretical models. The identifica- 
tion of the frequencies could be further improved by photo- 
metric studies to confirm this. Ultimately the sequencing of 
n-values could provide us with direct information about the 
stellar interior. 



For the reasons above we propose this star as an excel- 
lent candidate with which to test asteroseismic models and 
potentially give us insight into the pulsational behaviour of 
7 Doradus stars. 
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(b) x 2 values for 1 = - 3 for f p2 = 1.1863d- 1 . 
Best fit value is for (l,m) = (1,1). 
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(c) x 2 values for 1 = - 4 for f p3 = 1.6812c!- 1 . 
Best fit value is for (l,m) = (1,1). 
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(d) x 2 values for 1 = - 3 for f p4 = 1.215yd" 1 . 
Best fit value is for (Z,m) = (1,1). 
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(c) x 2 values for 1 = - 3 for f p5 = 1.5596d- 1 . 
Best fit value is for (l,m) = (1,1). 

Figure 6. Lowest x 2 values for each possible (I, m) combination 
for the final identified frequencies. The best fit x 2 is identified in 
bold. 
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8 SUPPORTING INFORMATION 

Additional Supporting information may be found in the on- 
line version of this article: 

Data file. The line profile data including Modified Ju- 
lian Date, velocity (on a relative scale) and the intensity 
at each of the 120 velocity sampling points for each profile. 
Line 1 is the axis scale, lines 2-478 are MJUO observations, 
lines 479-526 are La Silla observations and lines 527-586 are 
McDonald observations. 
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